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Abstract 


The radiation chemistry, thermal stability, and vapor pressure of solid-phase 
carbonic acid (H2CO3) have been studied with mid-infrared spectroscopy. A new 
procedure for measuring this molecule's radiation stability has been used to obtain 
intrinsic IR band strengths and half-lives for radiolytic destruction. Results are compared 
to literature values. We report, for the first time, measurements of carbonic acid's vapor 
pressure and its heat of sublimation. We also report the first observation of a chemical 
reaction involving solid-phase carbonic acid. Possible applications of these findings are 
discussed, with an emphasis on the outer Solar System. 


Key Words: Ices, IR Spectroscopy; Satellites, Surfaces; Cosmic Rays 


3 



the closest planet. In each case, carbonic acid may be formed, and for Callisto a tentative 
detection of H 2 CO 3 already has been made (Johnson et ah, 2004). 

To assess the formation and stability of carbonic acid in the Solar System, it is 
important to investigate the molecule's physical and chemical properties, but little such 
work has been published to date. Gerakines et al. (2000) compared the yields of H 2 CO 3 
made by exposing H 2 O + CO 2 ice mixtures to ion irradiation (~ 1 MeV H + ) and to UV 
photons (~ 10 eV). The same researchers measured carbonic acid’s intrinsic IR band 
strengths by the growth of products resulting from UV destruction of H 2 CO 3 . Earlier 
work also showed qualitatively that carbonic acid's vapor pressure is lower than that of 
H 2 O, CO 2 , and the observed reaction products, since H 2 CO 3 is the last of these to sublime 
under vacuum in the 200 - 250 K region (Moore and Khanna, 1991). A white color is 
likely for H 2 CO 3 made by acid-base chemistry (photographs in Loerting et al., 2000), and 
the work by Winkel et al. (2007) showed that the x-ray powder pattern of frozen H 2 CO 3 
is featureless. To our knowledge, little else, if anything, is known of solid-phase carbonic 
acid from experiments. 

In this paper, we reinvestigate the intrinsic IR band strengths of H 2 CO 3 and, for 
the first time, measure this molecule's radiolytic destruction at several temperatures. 

These new radiation experiments take into account amorphization of the sample. 
Furthermore, the highest temperature at which destruction measurements are made has 
been raised from ~10 K to 200 K. Temperature-dependent changes in the position and 
width of the H 2 CO 3 feature at 2618 cm H (3.82 pm) have been recorded. The first 
measurements of the vapor pressure and heat of vaporization of pure H 2 CO 3 are given, 
along with the first example of a low-temperature acid-base reaction of the molecule. 
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Changes in the IR spectra of irradiated ices were followed by Fourier- transform 
infrared (FTIR) spectroscopy using a Nicolet Nexus 760 instrument. In this setup, the 
incident IR beam passed through the sample, was reflected by the underlying aluminum 
mirror, and then passed through the ice a second time, and to the IR detector, for what are 
sometimes called transmission-reflection-transmission spectra. Measurements were 
made at 2-cnT 1 resolution from 5000 to 650 cm -1 , averaged over 150 scans. 

For studying the vapor pressure of H 2 CO 3 , the compound first was made by an 
acid-base reaction between a 1 molar solution of HBr (Sigma- Aldrich) and a 0.1 molar 
solution of KHCO 3 (Sigma- Aldrich), similar to the technique of Hage et al. (1993). A 
few microliters of the KHCO 3 solution were injected through a septum, using a syringe, 
onto a KBr substrate at 1 0 K, attached to the tail section of a closed-cycle helium 
cryostat. Next, a few microliters of the HBr solution were injected the same way to form 
a layer atop the frozen KHCO 3 solution. This process was repeated about ten times to 
increase the ice's thickness. Subsequent wanning of the sample to ~ 200 K removed the 
H 2 O and initiated a reaction between HBr and KHCO 3 to form H 2 CO 3 , changes that were 
followed with IR spectroscopy. The sample then was heated to 240 - 255 K, and IR 
spectra recorded over time, with a focus on the 1300 and 1500 cmf 1 features of H 2 CO 3 . 
Band areas were measured and combined with intrinsic band strengths, so called A 
values, to determine the vapor pressures and enthalpy of sublimation (AH SU b) of H 2 CO 3 
(Khanna et al., 1990). These measurements were made with a Mattson Polaris 
spectrometer operating in a conventional transmission mode. 
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3. 2. Radiolytic destruction 

The destruction of crystalline H 2 CO 3 by 0.8 MeV protons was followed by 
measuring the decrease in IR band areas after various doses. As an example, Figs. 3a and 
3b compare spectra in the 2900 - 1500 cm ” 1 (3.45 - 6.67 pm) region before and after 
irradiation to a dose of 2.0 eV molecule” 1 . The H 2 CO 3 bands are seen to decrease, 
indicating a loss of molecules, and at the same time H 2 O and CO 2 are formed (H 2 O is not 
shown). In addition, irradiation caused the H 2 CO 3 bands to widen, indicating 
amorphization of the crystalline sample. 

To accurately quantify carbonic acid's radiolytic destruction it was necessary to 
distinguish between spectral changes caused by (a) loss of H 2 CO 3 molecules and (b) 
amorphization. This was done by warming the sample to 200 K after each irradiation 
step to fully recrystallize the partially-amorphous ice and to sublime away the H 2 O and 
CO 2 formed by radiolysis. The ice then was recooled to 14 K, as shown in Fig. 3c, for 
comparison to the original spectrum of the unirradiated ice. Fig. 3a. The spectra of Figs. 
3a and 3c are similar, but the latter has slightly smaller H 2 CO 3 bands, caused by the 
destruction of crystalline H 2 CO 3 . 

The normalized band areas for F 12 C 03 have been plotted in Fig. 4 as a function of 
radiation dose. Table 2 lists all bands that were averaged for this graph along with their 
integration limits. Also in Fig. 4 are linear regression lines through the data points. The 
corresponding half-life doses for H2CO3 irradiated at 14, 100, and 200 K are then 11, 11, 
and 7 eV molec” 1 , respectively. 
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temperature. After this annealing cycle, some CO 2 often remained trapped in the H 2 CO 3 
(Fig. 3c). In order to relate only the amount of CO 2 formed to the amount of H 2 CO 3 
destroyed, at each radiation step we subtracted the band area of any remaining CO 2 in the 
annealed ice from the area of the CO 2 band recorded after the next irradiation. Table 3 
gives the results of these /ftTFCCft) measurements at 14 and 100 K, corrected for 
amorphization. Because most of the CO 2 product immediately sublimed away upon 
formation at 200 K, no band strengths were determined at that temperature. Note that the 

1 n I 

measurements in Table 3 are based on A( CO 2 ) = 7.6 x 10 cm molec (Gerakines et al., 
1995), and that no decomposition of H 2 CO 3 into CO appeared to occur. A few 
experiments with H 2 13 C 03 were conducted to verify that all of the CO 2 formation 
observed in our work was due to the proton irradiation, and not from leaks in the vacuum 
system. No such contamination was detected in any experiment. 

3.4. Radiation yield of H 2 CO 3 

The radiation-chemical yield, denoted G, of a substance is the number of 
molecules produced by absorption of 100 eV. Previously-reported values for G^COs) 
from H2O + CO2 (1:1) ices at 14 K were 0.028, 0.030, and 0.02 for MeV protons and UV 
photons (Gerakines et al., 2000), and for 10 keV electrons (Hand et al., 2007), 
respectively. These values were based on the growth of H2CO3 IR bands as a function of 
radiation dose, and represent the formation of H2CO3 within an amorphous ice mixture 
dominated by H 2 O and CO 2 . We repeated this type of experiment by irradiating H 2 O + 
CO2 (1:1) at 14 K and 50 K in small steps, the 1500 cm "" 1 band's area being measured 
after each irradiation. The column density of H2CO3 was calculated from equation (1) 
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and this may account for the 8 - 18 cm -1 shift of some bands with respect to their 

positions in the radiation-formed H2CO3. 

Figure 6 shows the decrease in the normalized average areas for the 1300 and 

1 500 cm'” 1 bands of H2CO3 as a function of time at five different temperatures. Each 

2 

decrease is related to a change in column density, the number of molecules per cm 
leaving the ice surface as a function of time. To determine column densities we used A 
values measured at 1 00 K. The slopes of similar non-normalized plots gave the 
sublimation fluxes (molec nf 2 sec ’) at each temperature. Equation ( 4 ) then was used to 
calculate the vapor pressure, p. 

Sublimation Flux = p / ( 27 imkT ) 1/2 ( 4 ) 

In ( 4 ), m is the mass of an H2CO3 molecule, k is the Boltzmann constant, and T is the 

'y 

absolute temperature, giving a vapor pressure in Newtons m , which was converted to 
units of bar. A plot of the calculated vapor pressure from 238 to 256 K is shown in Fig. 
7 a. Figure 7 b shows the same data graphed as ln(p) versus 1 /T, from which the slope 
gives the heat of sublimation as AH su b = 65.2 kJ mol -1 . 

3 . 6 . Chemical destruction 

In addition to measurements of both the sublimation and the radiolytic destruction 
of solid FI2CO3, we also have observed FI2CO3 loss by chemical reaction. Previously, we 
found that ammonia (NH3) blocks H2CO3 formation in irradiated solid-phase H2O + CO2 
+ NFI3 mixtures (Gerakines et al., 2000 ). In separate experiments, we have irradiated 
layered samples consisting of a mixture of solid H2O + CO2 over a layer of NH3, both 
ices being formed at ~ 10 K. Subsequent irradiation produced EI2CO3 in the upper layer, 
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5 shows the results of such a calculation, giving a half-life for carbonic acid at 100 K on 
both Europa and Callisto. 

The spectra we have recorded, such as in Fig. 1, illustrate the differences between 
carbonic acid in an amorphous matrix and pure crystalline H2CO3. On warming from 14 
to 240 K (Figs, lb and lc), some peaks shift, some bands narrow, and some splittings are 
observed. As an example, the broad, weak band near 2555 cm " 1 (3.914 pm) sharpens 
considerably and moves to 2612 cm " 1 (3.828 pm) on warming to 240 K. Figure 2 shows 
that this same feature then displays small, reversible shifts in position as the temperature 
of the carbonic acid is changed. The importance of documenting such spectral variations 
is demonstrated by Fig. 9, which overlays this same OH stretching feature of carbonic 
acid on reflectance spectra of Callisto and Europa. The spectrum of pure crystalline 
H2CO3 is shown as is one in which carbonic acid is trapped in an amorphous ice mixture. 
It is seen that the band shape and position for the amorphous ice provides the better fit. 
Any similar feature on Europa is within the noise of the data as demonstrated in Fig. 9. 
See Johnson et al. (2004) for a suggestion of H2CO3 as a possible candidate for Callisto. 
Additional details on the shifts and intensity changes of H 2 CO 3 features can be found in 
Winkel et al. (2007). 

Our analysis of data from irradiated H 2 CO 3 considers both radiation-induced 
chemistry and radiation-induced amorphization. Separating these effects is important 
because measurements of both H 2 CO 3 loss and CO 2 growth are needed for an accurate 
determination of intrinsic IR band strengths of carbonic acid. Along these lines, the band 
strengths we report in Table 3 are significantly different (> 50%) from some of the older, 
uncorrected values. 
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(25.2 kJ moF 1 ), which will assist H 2 CO 3 in remaining on planetary surfaces after the 
sublimation of the other two molecules. For comparison to other carboxylic acids, AH SU b 
is 62.5 kJ moF 1 for formic acid (HCOOH) and 67.9 kJ moF 1 for acetic acid 
(CH 3 COOH). See Calis-Van Ginkel et al. (1978). 

Having discussed our results, it is appropriate to point out some limitations and 
possible future work. Two sources of error in Table 1 are the unknown density and 
radiation stopping power of a 1 : 1 H 2 0 + C0 2 ice. Our approach was simply to assume 
these quantities to be an average of the values of the individual components. A direct 
measurement, particularly of the density, is desirable. 

The H2CO3 formation we report is for this molecule generated in an amorphous 
mixture of H 2 0 + C 0 2 ( 1 : 1 ). However, to quantify carbonic acid production we were 
forced to use our A values for crystalline H2CO3. The reason for this is that neither of the 
synthetic techniques we used to prepare H2CO3 resulted in the pure amorphous material. 
To our knowledge, pure amorphous H2CO3 has not yet been made and so no spectra or 
band strengths are available. A related difficulty concerned the purity of the carbonic 
acid in our vapor pressure measurements. The acid-base reaction used to make H2CO3 
gave KBr as a by-product. This is not expected to alter the vapor pressures of Fig. 7, but 
a check is desirable, such as with H2CO3 made by a different method. 

We also note that our vapor-pressure work was done with an IR spectrometer 
operating in a conventional transmission mode, while measurements of radiolytic 
destruction utilized reflection spectroscopy. The spectra in the two cases were essentially 
identical, as seen in Fig. 5. 
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Figure 5. A comparison of the IR spectra of H 2 CO 3 formed by the irradiation of an H 2 O 
+ CO 2 (1:1) ice and the acid-base reaction of warmed HBr + KHCO 3 . The upper 
spectrum was taken with the reflection method described in the text, while the lower 
spectrum was recorded in a conventional transmission mode. 

Figure 6 . Normalized areas (averaged) of the 1300 and 1500 cm '" 1 bands of H 2 CO 3 
plotted as a function of time for five different temperatures. Each point is from the 
average area of the IR bands at 1300 and 1500 cm -1 . 

Figure 7. The vapor pressures of H 2 CO 3 at five temperatures are shown. The slope of 
the regression line in the lower graph gives the heat of sublimation, AH s „b = 65.2 kJ 
mol -1 . 

Figure 8 . The IR spectrum of the H 2 CO 3 + NH 3 solid-phase reaction product compared 
to room-temperature ammonium carbonate, (NFLO 2 CO 3 . 

Figure 9: The IR reflection spectra of Europa and Callisto compared to spectra of 
crystalline H 2 CO 3 at 140 K and H 2 CO 3 mixed with H 2 O and CO 2 at 150 K. The H 2 CO 3 
absorbance spectra were inverted and arbitrarily scaled for this comparison. Callisto’s 
observed 3.880 pm (2577 cnT 1 ) feature is more similar to that of EI 2 CO 3 held in the 
amorphous ice. 
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Table 1 


Physical properties of ices 


Ice 

Molecular Mass 
(g mol' 1 ) 

Density 
(g cm" 3 ) 

Proton Stopping Power 
(MeV cm 2 g"') c 

h 2 o 

18 

1 

273 

C0 2 

44 

1.7 

240 

H 2 0 + C0 2 (1:1) 

31 

1.35 a 

256. 5 a 

h 2 co 3 

62 

l b 

254 


a Average value for a H 2 O + CO 2 (1:1) mixture. 
b Assumed value. 

Calculated for 0.8 MeV protons, according to method of Zeigler et al. (1985). 
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Table 3 


IR band positions and strengths (A) for H 2 CO 3 at 14 and 100 K 


Band Position 


^( 10' 17 

cm molec 1 

) 

Wavenumber 

Wavelength 





(cm^ 1 ) 

(pm) 

14 K 

100 K 

18 K a 

185 K b 

2749 + 2833 

3.634 + 3.530 

5.3 

3.3 

9.8 

— 

2618 

3.820 

7.5 

7.2 

16.0 

— 

1695 

5.900 

10.8 

14.8 

11 

35 

1503 

6.653 

5.2 

9.1 

6.5 

11 

1303 

7.675 

8.5 

12.3 

10 

12.6 

1038 

9.634 

0.15 

0.17 

0.14 

— 

875 

11.43 

4.5 

2.8 

5.6 

— 


a From 18 K photodissociation of H 2 CO 3 (Gerakines et al., 2000). 

b From 185 K implantation of H + into CO 2 ice to form H 2 CO 3 (Garozzo et ah, 2008). 


37 



Table 5 


Half-lives for crystalline H 2 CO 3 , based on the 100 K destruction rate in the laboratory, 
corrected for amorphization 


Environment 

Depth 

(pm) 

Volume Dose 
Rate 

(eV molec~' s -1 ) 

Half-Life 

Laboratory 

1.0 

1.3 x 10 ~ 3 

2.4 h 

Europa 8 

100 

1.0 x 10~ 8 

35 y 

Callisto 3 

100 

2.5 x 10 " 11 

1.4 x 10 4 y 


a Volume dose rates for Europa and Callisto from Cooper et al. (2001), 
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